perm filename TABLE.SAI[1,VDS] blob sn#498112 filedate 1980-02-03 generic text, type C, neo UTF8
COMMENT ⊗   VALID 00002 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002	BEGIN "TABLE"
C00006 ENDMK
C⊗;
BEGIN "TABLE"

	REAL ARRAY A[1:3,1:6],B[1:3,1:6],R[1:3],T[1:3,1:3],L[1:3,1:6],S[1:6];
	REAL LI,A1,B1,G1,X1,Y1,Z1,DX,DY,DZ,ALTHA,BETA,GAMMA,PI;
	DEFINE CRLF="'12&'15";

	PI←3.14159265
	OUTSTR ("TYPE IN A1,B1,G1 ");
       INCHWL A1,B1,G1;
	READ (5,205)X1,Y1,Z1
	WRITE (6,220)
	WRITE (6,230) A1,B1,G1
	WRITE (6,235)
	WRITE (6,240) X1,Y1,Z1
	DO 5 J=1,6
	READ (5,205) (A(I,J),I=1,3)
5	CONTINUE
	DO 10 J=1,6
	READ (5,205) (B(I,J),I=1,3)
10	CONTINUE
	ALPHA=0
	BETA=0
	GAMMA=0
	READ (5,205) DX,DY,DZ
	READ (5,200) ALPHA,BETA,GAMMA
	WRITE (6,245)
	WRITE (6,250) DX,DY,DZ
	WRITE (6,255)
	WRITE (6,260) ALPHA, BETA,GAMMA
	R(1)=DX+X1
	R(2)=DY+Y1
	R(3)=DZ+Z1
	ALPHA=(ALPHA+A1)*(2.*PI)/360.
	BETA=(BETA+B1)*(2.*PI)/360.
	GAMMA=(GAMMA+G1)*(2.*PI)/360.
	C1=COS(ALPHA)
	C2=COS(BETA)
	C3=COS(GAMMA)
	S1=SIN(ALPHA)
	S2=SIN(BETA)
	S3=SIN(GAMMA)
	T(1,1)=C3*C2
	T(1,2)=(C3*S2*S1)-(S3*C1)
	T(1,3)=(C3*S2*C1)+(S3*S1)
	T(2,1)=S3*C2
	T(2,2)=(S3*S2*S1)+(C3*C1)
	T(2,3)=(S3*S2*C1)-(C3*S1)
	T(3,1)=-S2
	T(3,2)=C2*S1
	T(3,3)=C2*C1
	DO 60 J=1,6
	DO 50 I=1,3
	L(I,J)=(T(I,1)*A(1,J)+T(I,2)*A(2,J)+T(I,3)*A(3,J))+R(I)-B(I,J)
50	CONTINUE
60	CONTINUE
	DO 70 J=1,6
	S(J)=SQRT(L(1,J)*L(1,J)+L(2,J)*L(2,J)+L(3,J)*L(3,J))
70	CONTINUE
	WRITE (6,265)
	WRITE (6,270)
	DO 80 J=1,6
	WRITE (6,275) J,S(J)
80	CONTINUE
200	FORMAT (3F6.3)
205	FORMAT (3F5.3)
220	FORMAT (' ','THE COMPONENTS OF INITIAL ANG. POSITION IN DEG. ARE')
230	FORMAT (' ','ALPHA=',1F7.3,2X'BETA=',1F7.3,2X,'GAMMA=',1F7.3)
235	FORMAT ('0','COMPS OF THE INITIAL POS. OF MOV FRAME ORIGIN ARE')
240	FORMAT (' ','X=',1F6.3,2X,'Y=',1F6.3,2X,'Z=',1F6.3)
245	FORMAT ('0','FOR DISPLACEMENTS FROM THE INITIAL POS. OF;')
250	FORMAT (' ','X=',1F6.3,2X,'Y=',1F6.3,2X,'Z=',1F6.3)
255	FORMAT ('0','AND COMPONENT VECTORS OF;')
260	FORMAT (' ','ALPHA=',1F7.3,2X,'BETA=',1F7.3,2X,'GAMMA=',1F7.3)
265	FORMAT ('0','THE LENGHTS OF THE STRAIN ELEMENTS WILL BE;')
270	FORMAT ('0',T30,'ELEMENT NUMBER',10X,'LENGTH,IN.')
275	FORMAT (' ',T38,I1,18X,1F6.3)
	STOP
	END